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We discuss a possibility of observing superfluid phenomena in a quasi-lD weakly interacting Bose 
gas at finite temperatures. The weakness of interaction in combination with generic properties of 
ID liquids can result in a situation when relaxational time of supercurrent is essentially larger than 
the time of experimental observation, and the behavior of the system is indistinguishable from that 
of a genuine superfluid. 
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Q^ . The discovery of Bose-Einstein condensation in ultra-cold dilute gases has opened up a unique possibility for the 
On study of fine aspects of superfluidity, including the specifics of lower dimensionality. Primarily, the new opportunities 
, come from the availability of isolated systems, free of contact with walls or any sort of substrate. A special interest is 
bJ[)' associated with the prospects of creating one-dimensional systems. Actually, one deals here with a quasi-lD situation, 
^ ' when the configuration of a trap yields toroidal spatial particle distribution with small transverse localization size 
' l±. At low temperature, and small interaction energy per particle, eo, as compared to the transverse- motion level 
] spacing, the system kinematics and collective dynamics are purely one-dimensional, though the value of the effective 
psj interaction is determined by 3D pair collisions (see below). 

] Strictly speaking, the ID superfluidity cannot exist at T ^ 0. - At flnite temperature the relaxational time for 
supercurrent is finite, i.e., independent of the system size in the limit L oo. Nevertheless, as will be shown below, 
analyzing the situation with a dilute gas one finds that the relaxation rate of the supercurrent is exponentially reduced 
■ with decreasing interaction, and that in a quasi-lT) weakly interacting system this can lead to a new quality when at 
] a finite temperature the lifetime of a supercurrent state considerably exceeds the lifetime of the isolated gas sample. 
OO ' We start with rendering basic notions and results for the low-energy properties of ID (super)fluids. First, we 
] remind that in ID there exists a condition on superfluid density, n^, and compressibility, k, determining whether the 
superfluid groundstate is stable against either single impurity, or disorder, or an external commensurate potential. 
This condition requires that the dimensionless parameter (m is the particle mass, h—1) 
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: K^l,^^^ (1) 

^ ; 

Ch , be less than corresponding critical value Kc (here c = ^JnJJmK is the sound velocity). For the most important for 

^ ■ our purposes cases of single impurity and random potential the critical values of K are K^. = 1 Q and — 2/3, 

• • , respectively Q. As we will see below, the case of quasi-lD weakly interacting Bose gas corresponds to K <C 1, so 

_ ^ ■ that the stability of the superfluid groundstate is guaranteed. 

, The low-energy spectrum of any one-component ID system in the superfluid phase is exhausted by bosonic sound- 

■ like excitations and a collective state supporting non-zero supercurrent, J/ = lunsl jniL, labeled by an integer index 
/ (we consider a system of ring geometry) 

B^^Y. + E ^^4^^/ . (2) 

where 6j creates a boson with energy LOq = cq and momentum q — 2T:k/ L, where k is integer. We have used a standard 
short-hand notation to specify the current state - obviously d\di is either 1 or and d\di = 1. The energy of the 
current state is given by 

Ei = . 3 

mL 

Note, that although the energy of the current state scales like ~ the momentum associated with it is system-size 
independent, Pj — iimsl . 

Within the long- wavelength effective Hamiltonian [describing the superffuid state in terms of the phase fleld], 
H = J dx[{ns/2){d<^/dxY + ('«/2)<i>^], where the original Bose fleld is written as 5" = y^e**, bosonic excitations and 
current states (sometimes called "zero mode" states |Q]) correspond to the following decomposition 
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$ = + 2t:Ix/L , <j>[d(l)/dx)dx = , 

where a; is a coordinate along the circle. Thus / is nothing but a phase field circulation 

I ^ ^ <j>{d^/dx)dx . (4) 

The very notion of the superfluid phase means that / is a well-defined quantum number even when rotational 
invariance is violated. In higher dimensions at T < Tc relaxation times r/ are extremely long with exponential 
dependence on the system size, and for all relevant experimental time-scales, Tcxp, the condition Texp t/ is satisfied. 
Special arrangements are to be made to shorten r/ (e.g., weak-link or tunnel-j unction systems). To find relaxation 
times in ID we may use the transition Hamiltonian which gives explicitly transition amplitudes between states 
|/, {Nq}) and |/', {N'^}) for arbitrary rotational invariance breaking terms. For the case of a local perturbation V{x) 
with spatial dimension smaller than the correlation radius 1 / mc, which without loss of generality may be replaced by 

V{x) VoS{x-xq) , 

the transition Hamiltonian acquires a form 

ifint = ffFon^4d,,e*2.„(7^7').oA,,,(6) + h.c. . (5) 

Polaronic-type exponential operator A (6) in this expression equals to 

A//K&)=exp|^^:^^^(^|^y^'sgn(g)(6,-6t)e'^^-''| . (6) 

The exact value of the numeric coefhcient g in Eq. (||) is outside the scope of the long-wavelength treatment used to 
derive Hint (it is of order unity for weakly-interacting Fermi gas). 

Formally, the problem has reduced to the problem of particle dynamics on a ID lattice with Ohmic coupling to the 
oscillator bath environment (see, e.g., review ||^) after polaronic on-site transformation. The dimensionless Ohmic 
coupling parameter 

a = (/ - ir/K (7) 

is large in our case, a > 1, and "particle" /current dynamics is incoherent at any finite temperature. We may then 
immediately utilize the well-known expression for the transition probability between states / and /' which derives 
from the Golden-rule expression for Hint (see, e.g., Ref. |^]) 



,2^,^ 2y^n |r[i + a + »g,,V(2^r)]p ^^^^,/2T 

^ir + r[l + a]r[l/2 + a] 



where 



^ = 2aln(^), (9) 

n = 2TTaT , (10) 

e//' ^Ej- Ej, , (11) 

7 is a numeric coefficient of order unity, and eo plays the role of the high-energy cutoff (at higher energies the dispersion 
curve is no longer sound-like and physics is determined by single-particle processes). Equation (||) was first derived 
in Refs. |,§. 

To find the decay rate rf^ one has to sum Eq. (H) over /', but since we are interested in the parameter range a S> 1 
and r <C eo the dominant contribution comes from /' = /± 1. It follows then from Eq. (||) that ^ ^ An^nsI/mL <C T 
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even for Ej ~ T, and we may neglect current energy transfer to the bosonic environment in Eq. (||). Also, for a ^ 1 
we may write approximately r[l + a]/r[l/2 + a] « ^/a. The final result may be written then as 

TV" V 7^0 / 

We see that for large a the decay rate is severely suppressed at low temperatures T < Tq ^ €o/2tt. Physically, this 
effect reflects small overlap between different vacuum states corresponding to quantum numbers / and /'. In practice, 
for a 3> 1 there is a "kinetic crossover" between the superfluid and normal behavior, i.e., between frozen and fast 
current relaxation, at T ^ Tq. 

To estimate realistic parameters for ultra-cold atomic gases we have to start from the original quasi-lD geometry. 
Let the transverse motion is confined by the symmetric parabolic potential with frequency uj±. Assuming that the 
3D scattering length a is small as compared with the oscillator length l±_ = 1 / ^miv± , we may derive the potential 
energy per particle from the 3D relation 

Cq = / dr± (pQ{r±) = nlloS , (13) 

m J 

where tpo is the wavefunction of the transverse motion (here n is the ID particle density) . For the parabolic potential 
we have 

eo — n2auj± . (14) 
Purely ID kinematics requires two conditions to be satisfied 

T < w± , eo<^uj± . (15) 



The second condition requires that a specific gas parameter, na, be small: na 1. When deriving Eq. ( |12[ ) we 
assumed that T <C eo to guarantee that transition rates between the supercurrent states are negligible. Obviously, 
this is a more severe restriction on the temperature range than the first inequality in Eq. (p^. 

We now turn to the parameter K for a translationally invariant system. (Weak disorder does not change the value of 
K dramatically.) For the weakly interacting gas compressibility is given simply by inverse UeS, and, in the low-energy 
limit, Hs ^ n. Thus, 

Experimentally it causes no difficulty to choose nl± > 1, and to make index K very small or parameter a — 1/ K ^ \. 

Let us consider an example of sodium gas assuming ID density n = 10^ cm~^. Then, na k, 0.25 and eo ~ 0.5 x u}^_. 
For the magnetic trap frequency u)±_ = 3 x 10^ s^^, we find nl± « 30, and extremely small K < 0.01 (!). The crossover 
temperature for the frozen current dynamics is roughly To = eo/27r « 2.5 x 10~^ K, and within the experimental 
range. 

We note that the effect of supercurrent stability considered above is essentially collective in nature and non- 
perturbative. Our transition Hamiltonian is valid only if initial and final states are given in terms of sound-like 
bosonic excitations - its validity may be questioned if the energy change in the transition, ^, exceeds cq. This imposes 
the restriction that interaction is not macroscopically small, Ucs S> 27r^/m_L. Also, even formally, for ^ > eg we find 
very fast relaxation rates from Eq. (||) (for estimates one may roughly replace in Eq. ( |l2|) 2ttT — > ^). We would like 
thus to verify that with all the above parameters we still may neglect energy (and momentum) transfer between the 
current and bosonic modes. Since £^ — Ej — Ei-i = Air^n{I — 1/2) /mL, for L = 1 cm we find ^ sa 8 x I0~^(/— 1/2) K, 
that is in the temperature range between 4 x 10^^ K and To we still have ^ < T for thermal values of / and system- 
size dependence is not crucial. On another hand, even if the energy has been dissipated into a single bosonic mode 
(in fact roughly a 3> 1 bosons are emitted/absorbed in the transition), the momentum of the boson would be only 
q = 47r(/ — 1/2)/ KL and much smaller than the momentum change of the current state 27rn. 

Thus, a quasi-lD interacting Bose gas with realistic parameters can support a supercurrent at T 0. The stability 
of such a state is guaranteed by extremely large relaxational time as compared to any reasonable experimental time. 
In connection with the question of preparation of supercurrent state, it is worthwhile to note that for an equilibrium 
ID system at finite temperature there is a finite probability to be found in a supercurrent state. This naturally follows 
from the fact that the energy of a supercurrent state scales like l/L, so that at finite T the statistical ensemble involves 
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a large set of different numbers /. After a deep cooling, the system will be found in a state with some particular /, 
randomly varying from one experiment to another | |To[ |. 
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